Have we found conclusive evidence for dark matter through direct 

detection experiments? 
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Abstract: We employ a SUSY-model-independent method to examine the remaining evidence for the low mass dark 
matter. Using the XENONlOO's recent result of 224.6 live days x 34kg exposure and PICASSO's result published in 
2009, we have obtained a constrain of couplings, \a n \ < 0.6 and \a p \ < 1.0, corresponding to the spin-dependent cross 



section of a n < 5.6 x 10 cm and cr„ < 1.6 x 10 



Spin-independent isospin violating dark matter model also 



failed to reconcile the recent result from XENON100 with the positive results from DAMA and CoGeNT. 

Key words: dark matter, spin-dependent, IVDM, XENON100, direct detection 
PACS: 95.35.+d, 14.80.Ly, 29.40.-n, 95.55.Vj 
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1 Introduction 

Since the first time the term "dark matter" was 
proposed by Fritz Zwicky in 1933Q, a variety of astro- 
physical and cosmological observations have provided 
convincing evidence, indicating that something invisible, 
however, having great influence on gravity does exist in 
our universe [2t5|- Now, it is generally accepted by astro- 
physicist that the ordinary atomic matter only makes a 
contribution of 5% to the universe, while the dark mat- 
ter (DM) takes the percentage of 23% and the remaining 
72% is dark energy [g. Meanwhile, a well-motivated 
extension of Standard Model, supersymmetric (SUSY) 
theory, provided excellent candidate for DM particles. 
It's now the most popular Weakly Interacting Massive 
Particles (WIMPs) among other non-standard WIMP 
models. 

In light of WIMP's property that it has weakly 
or less-than-weakly interaction with ordinary matter 
and the precondition, supported as well by a recent 
observation 0], that the Milky Way is embraced in the 
dark halo, DM can be searched for directly on earth. 
The terrestrial experiments designed for DM direct de- 
tections are all based on these preconditions. The energy 
deposited on the target after each hit by a passing DM 
particle will be transferred to detectable signals such as 
ionization, scintillation and phonon. 

The two modes of WIMP-nucleon interaction, which 
the direct detectors are aiming at, are spin-independent 
(SI) and spin-dependent (SD). SI coupling, the scalar 
mode, describes coherent interactions of the entire nu- 



cleus with WIMP. SD coupling, the axial one, describes 
the interaction depending on spin-content of the nucleus. 
Even though one of these experiments, DAMA/LIBRA, 
has been claiming success in finding signals of DM for 
more than a decade [H, 0], other experiments only give 
out null result and thus have made exclusion in the cross 



section and mass space for SI interaction[12j, [27J. Recent 
reports from CoGeNT and CRESST shows as well an 
evidence 13, 12 1- A survey flij] presenting consistence for 



these experiments in low mass area (10-20 GeV/c 2 ) in 
the SI mode drew a great deal of attention. The result 
of the survey applying the isospin- violating DM (IVDM) 
-nucleus interaction model, which assumes a different in- 
teraction strength between proton and neutron, however, 
should be checked by XENONlOO's updated resultsQ. 

In this article, a method of the calculation for SI 
WIMP-nucleon scattering will be briefly introduced in 
Sec. 2, which will be done by focusing on the IVDM 
model. In Sec. 3 we focus on the SD interaction to show 
the up-to-date constrains for coupling constants in the 
WIMP mass of 10 GeV/c 2 by analysing several leading 
experiments. Following that, we present the discussion 
and conclusion in Sec. 4. 

2 Spin-independent interaction 

For clastic WIMP-nucleon scattering, the event rate 
of an Earth-bound detector can be written as [13] 



R = N T 



Pd 
Mn 



dE h 



da 



vf(lT,vt)d 3 v (1) 
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where N T is the number of target nuclei, the local 
WIMPs density, p D =0.3 GeV/cm 3 ^, divided by the 
WIMP mass, M D , denotes the local number density of 
WIMPs. The integral interval of recoil energy E R are 
determined by experimental considerations. The lower 
limit of the velocity integral, u min = ^/ M A E R /2fi A , is 
the minimal value for a WIMP particle to deposit en- 
ergy E R on to a target atom, whose mass is M A , and 
/i_4 is the WIMP-nuclei reduced mass. The upper limit 
v eS c = 544 ± 64 km/s[18( is the local Galactic escape 
velocity. Maxwellian distribution, with characteristic 
velocity v = 220 km/s|18], for f(v) is assumed for a 
staple dark halo model[l3, where "if is the DM velocity 
onto the detector while v E is the earth velocity relative 
to the static galaxy. 

The differential cross section, da/dE R , is model- 
dependent, which can be written in a general form, 



da 
dE~ R 



(2) 



where the zero-momentum-transfer cross section a A and 
the form factor F 2 (q 2 ) are different for the two interac- 
tion modes. 

In SI mode[20|, the form factor proposed by Helm[2l| 
is applicable to various sorts or targets, and the general 
form of cross section can be 



= -^ A [f p Z + f n (A-Z)f 

IT 



(3) 



where f p and f n are the coupling constants for WIMP- 
proton and -neutron scattering, which for Majorana 
WIMPs is usually taken for equality, f p ~ f n = /. Thus 
Eq. (3) turns to erf 7 = (4/ 2 /-k)h 2 a A 2 , that is proportion 
to the square of nucleon number A. Otherwise, f p ^ f n 
is called the isospin- violating Dark Matter (IVDM). In 
Ref . (20| . the assumption of a rate, / p //„ = —0.7, suc- 
ceeded to break away from the constrain of XENON100 
and reconciled CoGeNT with DAMA. Here, we apply 
the latest result of XENON100 with the same IVDM 
coupling to make a exclusion curve and check the agree- 
able region of CoGeNT with DAMA in Ref. [20], find- 
ing that the once successful model failed the exam this 
time, shown in Fig.l. As XENON100 uses xenon as tar- 
get, which has seven isotopes, the set of f p /f„ = —0.7 
has reached its limit to reduce the constrains from the 
detector. It will not be able to invalidate XENONlOO's 
boundary any longer. 
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Fig. 1. Favored region by CoGeNT/DAMA and 
the exclusive curve by XENON100 for IVDM with 
/,//„ = -0.7. 



However, the set of f p /f n = —0.7 is specific. A more 
general method of deciding the cross section (the cou- 
pling constants) for WIMP's SI scattering provides an 
overall picture for a given WIMP mass, which is shown 
in Fig. 2. The up-to-date f n — f p region is constrained by 
XENONlOO's recent result as an orange ellipse. 
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Fig. 2. Spin-independent couplings for a 10 
GeV/c 2 WIMP 



In Fig. 2 the allowed coupling constants, /„ and /„, 
is the inner region of the ellipses. Colored bands rep- 
resent the allowed parameters from the CoGeNT and 



DAMA evidence. Although XENON100(2011)'s results, 
the green ellipse, left some room for CoGeNT and 
DAMA's evidence, shown in the purple regions, the re- 
cent XENONlOO's results, the orange one, has com- 
pletely excluded them. 



the total nuclear angular momentum. For single nucleon 
whose spin and total angular momentum are of the same 
value, J = 5 = ±1/2, 



24 „ 2 2 2 

• —Lr F fi p a p7l 



(10) 



3 Spin-dependent interaction 

Besides the well-motivated Majorana fermion, a su- 
perposition of bino, wino and higgsino appears to be 
recently the most popular model for WIMP candidate. 
The Majorana WIMP's interaction can be both SI, intro- 
duced in Sec. 2, and SD, to be discussed in this Section. 
The gauginos-higgsinos mixed WIMP (or the so called 
the lightest SUSY particle, LSP)-nucleon interaction 
is relevant to its components. Since the precise com- 
position of the LSP is unknown, it is preferable to use 
a SUSY- model-independent method [22| to do the survey. 



The general form of SD cross section is given as [22 



'(?) = 



32G 2 Ff i 2 A 
2J+1 



(4) 



where G F is the Fermi constant, fi A is the nucleus- WIMP 
reduced matter. When normalized, S(q)/S(0) is the 
form factor for SD mode, which is the counterpart of 
F 2 (q 2 ) in SI mode. The expression of S(q) can be ex- 
panded to 

S(q) = [alS P p(q) + a p a n S pn (q) + a 2 n S nn (q)} (5) 



S pp = S 00 (q) + S 01 (q) + S 11 (q) 
S nn = Soo (q) — Soi (q) + Su (q) 
S pn = 2[S 00 (q)-S 11 (q)] 



(6) 



S i: j (i,j = p,n) is the spin structure, that is specific to 
each nucleus. The coupling constants are denoted as a p 
and a n . In the case of zero momentum transfer, 



5(0): 



^±1aV ( j + i) 



(7) 



where A is given in the form[24j, 



(N\{a p Sp + a n S n )\N[ 



N 



J 



N 



= j(a p (Sp)+a n (S n }) (8) 



From Eq. (|4j) (jTj) (JSj) , the zero momentum cross section for 
SD is derived 



32 

a A = — G F fi' A [a p (S p ) + a„ (S n 



,J+1 
J 



(9) 



where (S p ) and (S n ) are the expectation of spin for pro- 
ton and neutron by the odd-group nuclear model, J is 



Thus Eq. 
section 



TU| can be expressed by single nucleon cross 



"5 



:J+1 
J 



(11) 



Owing to the relation of a p and a n , it is impossible to 
work out the single nucleon cross section(or its upper 
limit) with WIMP. Although we can set a ration of aj o„, 
it becomes SUSY-model-dependent. Experiments (2H . I2H 
l3~fl | choose to set a v to zero when calculating the WIMP- 
neutron cross section, and vice versa. This is done by 
assuming a A = <t a p , where 



32 

7T 

32 



G\n\(a v (S p )) 



cr A — — G 2 F n 2 A (a n (S n )) 



. J + l 
J 

2 J+l 



(12) 



7T - J 

Despite the ambiguous physical meaning, the a A and <r A 
are two useful auxiliary mathematical expressions, with 
the aid of which, the SD cross section can be rewritten 

as 

a A = (,M±V^f (13) 

the sign of which is identical with the sign of (S p ) / (S n ). 
The comparison between Eq. (fTU|)and Eq. (fT2j) gives 



1 J 3f4 



n.p 
A 



(14) 



After combining Eq. @ (HJ) © , we find the integration of 
Eq. (1) brings out an elliptical curve in the a p -a n coor- 
dinate system, which can be expressed as 



N- 



Aa 2 +Ba p a n - 



-Ca 2 



(15) 



Eq. (|15l) is a conic section, which can be an ellipse, 
whose center is the origin of coordinates, if B 2 < AAC. 
Otherwise it can be two open curves: hyperbola when 
B 2 > 4AC or two parallel lines if B = ±2y/AC. If the 
target contains only single nucleon that is sensitive to 
SD interaction, where being sensitive for SD interaction 
means the nucleon should be odd-neutron or odd-proton 
to have non-zero spin, the Eq. (ITBl changes to 

(16) 



N' = (A'a p + C'a n ) 2 



which shows that the ellipse has degenerated into two 
parallel lines. If Eq. (fT2|) is used to calculate the pure 
neutron and pure proton cross section (or their upper 
limits), actually, it's taking the two elliptical points on 
the coordinate axes, (a n , 0) and (0, a p ), to obtain the a p 



and <7 n varying with the WIMP mass. 

The XENON100 experiment uses the targets with an 
isotopes abundance of 26.2% for 12g Xe and 21.8% for 
129 Xe, both of which are sensitive to WIMP-neutron SD 
scattering. The result shows that two candidate events 
have been observed in the energy range of 6.6-30.5 keV nr 
from the 224.6 live days x 34 kg exposure [lj]. In the 
analysis, a set of energy range of 6.6-43.3 keV nr (3-30PE) 
is used with the expected background of 1.0 ±0.2. Ap- 
plying the Feldman- Cousins procedureflij], we obtained, 
up to the present, the strongest limit of a„. In the cal- 
culation, we apply the spin structure from Menendez, 
Gazit and Schwenk's work (MGS for short) [28| and an- 
other work by Ressell and Dean (RD)[29| is taken for 
reference. 

Since evidences from CoGeNT and DAMA match- 
ing in the WIMP mass of 10 GeV / c 2 flH| have not been 
checked in the SD mode, we calculated at this mass point 
to find out the agreeable region in the a p -a„ area. This 
can be done by the method introduced above or use an 
alternative way 22] : if the nuclear cross section (for ev- 
idence found) or its upper limit (for null result) is re- 
ported, which is obtained in the procedure that account 
the total events' rate for a single isotope (Eq. ©-(13)), 
we can use Eq. (fl~7|) straightly, 



± 



lim(Ai) 
r \ t y 



lim(Ai) 



(17) 



where <J p ln ^- A ^ is the proton (neutron) experimental up- 
per limit for a certain nuclide, which for evidence found 
experiments is the certain cross section. 



In the procedure of dealing with CoGeNT's 
spectrum [16J, we adopt the most stringent surface event 
rejection to check the lower limit of its couplings (or 
cross section) with an 99% confidence level. Other cou- 
plings' region under a milder surface rejection will surly 
be excluded if the most stringent one can even not sur- 
vive. 




Fig. 3. Spin-dependent couplings allowed by 
XENON100 and PICASSO at a WIMP mass of 
10 GeV/c 2 (filled in orange). 



We take PICASSO's[3lJ result for a p constrain since 
its low threshold energy qualifies it to provide an exam 
in the low mass (10 GeV/c 2 ) region. The experiment use 
the superheated C 4 F 10 to capture the droplet 's explosion 
into a vapor bubble after hitting by a WIMP. In virtue of 
its single nuclear interaction in SD mode, we can easily 
obtain the CT^ mA from the published er' lmA by 



(S n ) 



< imA (s P y 

The spin expectation is listed in Tabel[TJ 

Table 1. Spin values for relevant nuclides. 



(18) 



Nucleis 


Odd 


J 


(S P ) 


(Sn) 


Ref. 


19 F 


P 


1/2 


0.441 


-0.109 


[32] 


23 Na 


P 


3/2 


0.248 


0.020 


[29] 


73 Ge 


n 


9/2 


0.030 


0.378 


[33] 


129 Xe 


n 


1/2 


0.010 


0.329 


[28J 


131 Xe 


n 


3/2 


-0.009 


-0.272 


[28] 



In Fig. 3, the allowed region of a p and a n is con- 
strained both by XENONlOO's ellipse and PICASSO's 
dashed lines, which is filled in orange. CoGeNT and 
DAMA's evidence(colored bands) is excluded obviously. 

4 Discussion and conclusion 

In this article, we have examined the signal reported 
by DAMA/NAI and CoGeNT with the recent result from 
XENON100 and PICASSO experiments. The remained 
unchecked 10 GeV/c 2 DM evidence is excluded for both 
the SI IVDM model and the SD model. 

• The IVDM model with f p /f n = -0.7 is unable to 
keep the CoGeNT and DAMA's compatible region 
evading from the XENONlOO's constrain. That 
is to say, in SI mode, there's no evidence can be 
affirmed from DM direct detection experiments ei- 
ther. 



• With the constrains from XENON100 (2012) and 
PICASSO (2009), for 10 GeV/c 2 DM mass, we 
have obtained the allowed couplings' region of 
\a n \ < 0.6 and \a p \ < 1.0, corresponding to the cross 
section a n < 5.6xl0~ 38 cm 2 and a p < 1.6xl0~ 37 cm 2 . 
Thus, in SD mode, no WIMP signal is compatible 
within all experimental results. 

For SD calculation, a model-independent method has 
been introduced. While experiments are used to report 
the cross section in the pure neutron or pure proton 
way, it is insufficient to give a full comparison between 
different detectors. This has been shown in Sec. 3 for a 
10 Gev/c 2 WIMP. Other masses are also excluded using 
the same method. 

I'd like to thank my supervisor K.Ni for useful dis- 
cussion and helpful advices given on the improvement of 
this work. 
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